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ABSTRACT 

We have obtained / band Tully-Fisher (TF) measurements for 522 late-type galaxies in the 
fields of 52 rich Abell clusters distributed throughout the sky between ~50 and 200h~ 1 Mpc. 
Here we estimate corrections to the data for various forms of observational bias, most notably 
Malmquist and cluster population incompleteness bias. The bias-corrected data are applied to 
the construction of an / band TF template, resulting in a relation with a dispersion of 0.38 
magnitudes and a kinematical zero-point accurate to 0.02 magnitudes. This represents the most 
accurate TF template relation currently available. Individual cluster TF relations are referred to 
the average template relation to compute cluster peculiar motions. The line-of-sight dispersion 
in the peculiar motions is 341 ± 93 km s , in general agreement with that found for the cluster 
sample of Giovanelli and coworkers. 



Subject headings: galaxies: distances and redshifts — cosmology: observations; distance scale 



1. Introduction 



Deviations from smooth Hubble flow arise as a result of large scale density fluctuations. Qualitatively, 
a "convergence depth" is the distance out to which significant contributions to a galaxy's peculiar motion 
are made. In linear theory, the fraction of the peculiar velocity of a galaxy, contributed by the mass 
distribution within a radius R, can be written as 

V(R)^^1 J S(r)^W(R)d 3 r, (1) 

where the Hubble constant H Q = lOO/i km s _1 Mpc -1 , f2 is the cosmological matter density parameter, S 
is the mass overdensity, and W{R) is a window function of width R (a Gaussian or top-hat, for example) 
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centered at r=0 (Peebles 1993). If the distribution of matter approximates homogeneity on larger scales, 
then contributions to the peculiar velocity will eventually taper off with increasing distance. 

Observational estimations of the local convergence depth are facilitated through measurements of the 
reflex motion of the Local Group with respect to spherical shells of increasing radii. The first such program 
was carried out by Rubin et al. (1976) on an all-sky sample of 96 Sc galaxies between redshifts of 3500 
and 6500 km s _1 . They reported the motion of the Local Group with respect to the shell of galaxies was 
454±125 km s _1 towards (l,b) = (163°, —11°), significantly different from the apex of the CMB dipolc, 
suggesting a large bulk flow for the shell and a convergence depth significantly larger than ~50/i _1 Mpc. 
Subsequent work in this area has yielded conflicting results. 

The situation can be grossly expressed into two main views. A picture of a small convergence depth 
(cz < 5000 km s -1 ) was first suggested by Tammann & Sandage (1985) and successively by Dressier et 
al. (1987) and Lynden-Bell et al. (1988). This picture is also supported by early measurements of the 
distribution of IRAS and optical galaxies (Lahav, Lynden-Bell & Rowan-Robinson 1988; Lynden-Bell, 
Lahav & Burstein 1989; Strauss et al. 1992; Hudson 1993). It should be mentioned, however, that some 
IRAS dipoles suggest significant contributions to the Local Group motion arise from the distribution of 
objects between 5000 and 10,000 km s _1 (e.g. Rowan- Robinson et al. 1990). The latter studies profit from 
flux dipoles or a peculiar velocity field that is solved for iteratively, initially supposing that recessional 
velocities are indicative of distance and that light traces mass. The above observational programs that 
utilize this technique find 80% or more of the Local Group motion derives from matter within ~50/i _1 Mpc. 

Proposals for a large convergence depth (cz > 10,000 km s _1 ; Scaramella, Vettolani & Zamorani 1994; 
Tini-Brunozzi et al. 1995; Branchini & Plionis 1996; Plionis & Kolokotronis 1998), mainly resulting from a 
similar analysis of the distribution of clusters of galaxies, suggest that much of our motion may be produced 
by mass concentrations as far as ^13,000 km s _1 (e.g. the "Shaplcy Supercluster"). The latter view is 
corroborated by reports of bulk motion in the local Universe (within 6000 km s -1 ) via Tully-Fisher (TF) 
measurements by Willick (1990), Courtcau et al. (1993) and Mathewson, Ford & Buckhorn (1992), and by 
the recent analysis of the dipole of cluster brightest ellipticals by Lauer and Postman (1994, henceforth 
LP). The LP claim, based on the peculiar motions of all Abell clusters to 15,000 km s _1 , suggests the 
local volume of space within ^lOO/i -1 Mpc is traveling towards (l,b) = (343°, +52°) at 689 km s _1 . 
The disagreement between these competing views of the convergence depth scale is wide, with important 
implications for cosmological models, which are generally unable to accommodate bulk flows with scales as 
implied by the large convergence depth camp (e.g. Gramann et al. 1995). 

Fresh work from Hudson et al. (1999) and Willick (1999) challenge the direction of the LP bulk flow 
vector. Hudson et al. employ the fundamental plane relation for some 700 ellipticals in 56 clusters (3000 
£ cz £ 14,000 km s" 1 ) to find a bulk flow of 630±200 km s" 1 in the direction (l,b) = (260°, -1°). A 
different bulk motion of ^700±250 km s _1 , towards (I, b) = (272°, 20°) originating from TF measurements 
in 172 cluster galaxies (9000 <; cz <; 13,000 km s _1 ) and 72 other galaxies (cz < 30,000 km s^ 1 ), is claimed 
by Willick. Such large bulk flows are not consistent with other recent observational work. Riess, Press & 
Kirshner (1995) used 13 SN la observations and found evidence for a small local bulk flow. In a contribution 
based on TF distances of about 2000 galaxies within 9500 km s _1 , Giovanelli et al. (1998a) similarly report 
evidence against the existence of a large scale local flow, and support for a relatively small convergence 
depth. Many of the aforementioned studies cannot convincingly exclude the existence of large-scale bulk 
flows, either because their sampling is too sparse (Watkins & Feldman 1995) or of limited depth. This is 
particularly important in view of the claims that asymptotic convergence of the Local Group reflex motion 
may only be reached at distances well in excess of 10,000 km s^ 1 (Scaramella, Vettolani & Zamorani 1994, 



-3 - 



Tini Brunozzi et al. 1995, and Branchini & Plionis 1996). This work aims at the direct determination of the 
TF relation for an all-sky cluster set extending to distances exceeding the highest reported values of the 
convergence depth. 

In addition to ensuring that our TF relation is valid for usefully large distances, a second issue of 
concern regards the amplitude of systematic errors in the TF template. The Giovanclli et al. (1998a) 
sample of peculiar velocities obtained using the TF method is referred to a template relation based on 
the SCI sample, a collection of 782 TF measurements in the fields of 24 separate clusters between 1000 
and 9000 km s _1 (Giovanelli et al. 1997a, b; hereafter G97a,b). Though the relative proximity of the SCI 
allows a broad stretch of observable galactic properties and thus makes it an ideal sample to study several 
characteristics of the TF template relation, it has limits as to how accurately the relation's zero point can 
be pinned down. A larger and deeper cluster sample has two main advantages. First, the increased number 
of clusters reduces the impact of statistical "shot noise." In addition, since the magnitude offset produced 
in the TF diagram by a given peculiar velocity decreases with the target distance, the scatter produced by 
peculiar velocities of distant clusters about the template relation is reduced; thus they are better suited to 
determining the template relation's zero point. Here we present the results of a program designed to probe 
the large-scale peculiar velocity field to ^200/i _1 Mpc, consisting of spectroscopic and photometric data for 
an all-sky sample of 522 galaxies from the fields of 52 clusters between ~50 and 200h~ 1 Mpc (hereafter the 
'SCIF sample) .0 

The overall scatter of the I band TF relation of about one-third of a magnitude translates to an 
uncertainty of 15% in redshift-independent distance measurements. This means that for an individual 
galaxy at, say 15,000 km s _1 , the method is able to predict the distance to within 2250 km s _1 . This value 
is considerably larger than the typical value of peculiar velocities, which are of order 500 km s _1 or less. 
Even if 50 objects per cluster were to be measured, the distance of a cluster at 15,000 km s _1 would not be 
characterized to better than 2250/v / 50 ~ 300 km s _1 ; our measurements include many fewer galaxies per 
cluster field, typically 10, some of which turned out to be cluster members. The main purpose of this study 
was thus not the determination of accurate individual peculiar velocities of remote clusters, but rather to 
combine measurements in many clusters to obtain a global solution for the dipole of the velocity field. 

This is the fifth paper in our series on the local convergence depth. The observational data are 
presented in Dale et al. (1997, 1998, 1999; Papers I, II, and IV; the data can be obtained by contacting the 
first author). The core result of this work, that the local dipole flow to about 200h~ 1 Mpc is consistent 
with a null bulk motion, is described in Dale et al. (1999; Paper III). Details of the sample and its selection 
are covered in Section 2 of this work, while Section 3 covers the construction of the universal TF template 
relation. Results for the peculiar velocity sample are given in Section 4, and we summarize our findings in 
Section 5. 



2. Sample Selection 

Clusters of galaxies are used as increasingly pliant tools in cosmology. For our particular concern, the 
cluster peculiar velocity distribution reliably matches that of the underlying smoothed matter's velocity 



2 The SCI and SCII are complementary samples of cluster TF data. The SFI is a completely independent sample of TF data 
for some 2000 field galaxies. Details on the SCI and SFI samples can be found in Giovanelli et al. (1998a, b) and references 
therein. 
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distribution (Bahcall et al. 1994a, b; Gramann et at 1995). Furthermore, though observations of clusters of 
galaxies only sparsely sample the large scale velocity field, they can do so more accurately than individual 
galaxies can. This advantage arises because measurements of many galaxies within a single bound system 
can be made. As peculiar velocities are driven by gravitationally growing density fluctuations, precise 
comparisons of large scale peculiar velocities can be made with those predicted by cosmological theories 
(Watkins & Feldman 1995; Cen et al. 1994; Feldman & Watkins 1994; Strauss et al. 1995; Croft & 
Efstathiou 1994; Borgani et al. 1997). 



2.1. Sample Definition 

Clusters of galaxies are practical TF targets. Their outskirts are well populated by spiral galaxies, thus 
a small number of wide-field images with a modest-sized telescope can effectively "map" a cluster, and they 
will likely contain numerous TF candidates. 

A second consideration that suggests the use of clusters in TF experiments is the determination of the 



TF template. As will be outlined in Section 3.1, the template is preferentially measured within a cluster, 
or determined as an average of templates from measurements in separate clusters. The template slope is 
best measured if a broad dynamic range of TF parameters can be observed, preferentially obtained from a 
relatively nearby sample. In contrast, the kinematical zero point is better characterized at higher redshifts, 
where cosmic peculiar motions play a smaller role in shifting objects away from the template. Our motive 
of improving the TF zero point accuracy plays an important role in characterizing the redshift distribution 
of our sample. 



A third advantage to using clusters of galaxies in TF work is the ~vA r increase in statistical accuracy 
per cluster they afford if TV TF measurements per cluster are available - the estimate of a system's peculiar 
motion is more accurate when information from multiple objects is used, the distance differential between 
galaxies in a cluster being negligible for our sample. However, if an all-sky survey of peculiar motions is the 
goal, concentrating observations to a single cluster or to a small number of clusters parallels an increase in 
overall sample sparseness. The choice between a densely sampled volume of low accuracy peculiar velocities, 
and a sparse sample of accurate peculiar velocities ultimately depends on the particular goals of the study. 

We early on investigated what sample characteristics would facilitate the most accurate determination 
of the TF zero point and the local bulk flow. To ascertain the optimal distribution of clusters and number 
of galaxies per cluster to observe, we relied on numerical simulations. The only limitation enforced was the 
total number of objects that could be observed, for target accuracy and a reasonable project timescale. 
A wide range in the number of clusters and the number of objects per cluster was explored. A variety of 
models of the peculiar velocity field (multi-attractor, linear bulk flows, quiescent, etc.) was imposed, and 
different models for the shape of the Zone of Avoidance (ZoA) were considered. These studies profited from 
numerical simulations of various types of cold dark matter models kindly provide by S. Borgani (see Paper 
III for further details). TF errors were assigned using random Gaussian deviates of the scatter distribution 
in G97b. The results, which took into consideration that clusters of galaxies have positions that are 
correlated in space (and therefore a random choice of N clusters does not guarantee that N independent 
points are sampled), suggested that about 50 clusters of galaxies needed to be observed, with measurements 
of at least 7-8 TF galaxy distances in each. 
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2.2. Cluster Selection 

The initial characterization of the SCII cluster sample predates this study and originates from 
preparatory work done by R. Giovanclli for the SCI TF study. Clusters were selected using the Abell rich 
cluster catalog as a guide (Abell, Corwin & Olowin 1980; hereafter ACO). Expediency played an important 
role in choosing the parent sample of clusters, of which this sample is a subset. Nearby clusters (cz <, 
10,000 km s _1 ) with preexisting H I velocity widths and / band fluxes were favored, and the more distant 
clusters (cz ;> 10,000 km s _1 ) that already had a large number of redshifts available were likewise strongly 
considered as that availability improved the prospects for their kinematical characterization. Table [l] lists 
the main parameters of the 52 chosen clusters for this work. The parameters listed include: 
Column 1: Standard name according to Abell/ACO catalog. 

Columns 2 and 3: Adopted coordinates of the cluster center, for the epoch 1950; they are obtained from 
ACO, except for the entries A1983b and A2295b, systems found to be slightly offset from A1983 and A2295 
in both sky position and redshift. 

Columns 4 and 5: systemic velocities in the heliocentric and in the cosmic microwave background reference 
frame, respectively, where we assume the motion of the Sun with respect to the CMB is 369.5 km s _1 
towards (I, 6)=(264?4,48?4) (Kogut et al. 1993). For all the clusters we derive a new systemic velocity, 
combining the redshift measurements available in the NED^ database with our own measurements. An 
estimated error for the systemic velocity is parenthesized after the heliocentric figure. 
Column 6: the number of cluster member redshifts used in determining systemic velocities. 
Column 7: Abell richness class. 

Column 8: Bautz-Morgan code (Bautz & Morgan 1970) as listed in the ACO. 

Figure [l] displays the sample in Galactic coordinates. The symbol sizes are inversely proportional to 
the cluster redshifts; two examples are given in the lower left for scale. An alternative display of the sample 
is shown in Figure ^|, a stereographic view of the sample in Galactic Cartesian coordinates. The dashed 
circles have radii of 5, 10, and 15 thousand km s _1 and the solid lines in the [X,Z] and [Z,Y] plots are for 
|6| = 20° and identify the ZoA. This distribution recalls our first criterion in selecting the sample: The 
data set should uniformly sample as much of the sky as prudently feasible. Since the main thrusts of this 
work are to recover a bulk flow measurement and to accurately determine a kinematical offset even in the 
presence of such a flow, an all-sky sample is required. Unfortunately, the paucity of clusters and the large 
and uncertain Galactic extinction in the direction of the ZoA prohibit us from sampling that portion of the 
sky. In addition, the likelihood of foreground star contamination increases dramatically for objects at low 
Galactic latitudes. Our formal criterion was to select clusters at |6| J> 20° from the whole sky. 

A histogram of the SCII redshift distribution is presented in Figure ||. The range of redshifts runs from 
5000 km s" 1 to 25,000 km s" 1 , with fully 90% of the clusters lying between 7000 km s" 1 and 19,000 km s _1 . 
The average CMB redshift of the SCII clusters is 12,050 km s _1 when clusters are weighted according to 
the square root of the number of TF measurements available. It is evident from Figure |^ that our sample's 
distance range provides an opportunity to effectively test claims of bulk flow motions on scales of 100/i -1 
Mpc. 



3 The NASA/IPAC Extragalactic Database is operated by the Jet Propulsion Laboratory, California Institute of Technology, 
under contract with NASA. 
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Table 1: The SCII Cluster Sample 
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064) 
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Table @ (Continued) 



i^l- , ,-.J- 

(cluster 


K. A. 


Dec. 


(l,b) 


CZhclio 




AT 


n 


ts-M 




B1950 


B1950 


B1950 


km s _1 


km s _1 








A3667 
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A3716 
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(346,-39) 
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1 


I-II: 


A3744 
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11386(089) 
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42 


1 


II-III 


A2457 


223312 


+011300 


(069,-47) 


17643(110) 
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20 


1 


I-II: 


A2572 


231554 


+182800 


(094,-39) 


11857(100) 


11495 


38 





III 


A2589 


232130 


+163300 


(095,-41) 


12288(095) 


11925 


58 





I 


A2593 


232200 


+142200 


(093,-43) 


12414(086) 


12049 


67 





II 


A2657 


234218 


+085200 


(097,-50) 


12028(137) 


11662 


40 


1 


III 


A4038 


234506 


-282500 


(025,-76) 


9012(063) 


8713 


180 


2 


III 




Fig. 1. — The all-sky distribution of the SCII sample in Galactic coordinates. The circles are drawn inversely 
proportional to the cluster redshifts. The two examples in the lower left give the scale. Asterisks mark the 
apices of the motion of the Local Group with respect to the CMB and the LP cluster inertial frame. 

2.3. Galaxy Selection 

To determine locations of target fields to be imaged, we visually scanned the Palomar Observatory Sky 
Survey plates for regions in the clusters containing promising disk systems appropriate for TF work. The 
selection of target galaxies for this study stemmed directly from the images obtained. A discussion of the 
imaging for this project is contained in Papers I, II, and IV. Each image in clusters chosen for spectroscopy 
was searched for spiral disks with the following properties: 

(i) disk inclinations ;> 40°; 

(ii) lack of dominating bulges - bulgy disk systems tend to be gas deficient and thus undesirable for emission 
line spectroscopy. Moreover, morphological homogeneity is preferable, to limit the effects of morphological 
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Fig. 2. — The stereographic representation in Galactic Cartesian coordinates of the 52 clusters of galaxies 
that make up the SCII sample. The dashed circles have radii of 5000, 10000, and 15000 km s _1 . The solid 
lines in the [X,Z] and [Z,Y] plots indicate |6| = 20°, the approximate extent of the Zone of Avoidance. 



bias (Section 3.3.5) 



(iii) no apparent warps/interacting neighbors; and 

(iv) no nearby bright stars which may affect flux measurements. 

It should be noted that the above properties served as guidelines for selecting TF candidates, but 
occasionally the guidelines were not strictly followed, as occasionally our hand was forced by the vagaries of 
telescope allocations and weather conditions. Coordinates and position angles for approximately 2250 TF 
candidates were obtained from the Digitized Sky Survey^ and are accurate to better than 2". The full list 
is available from D.A.D. upon request. 



4 Thc Digitized Sky Surveys were produced at the Space Telescope Institute under U.S. Government grant NAG W-2166. 
The images of these surveys are based on photographic data obtained using the Oschin Schmidt Telescope on Palomar Mountain 
and the UK Schmidt Telescope. 
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Fig. 3. — The redshift distribution of the clusters. 



A TF display of the collective raw data can be seen in Figure ||. The term "raw" implies that the 
absolute magnitudes Mj and the rotational velocity widths W are corrected for all the effects described 
in Papers I and II (extinction, inclination, etc), but the effects of cluster incompleteness bias and peculiar 
motion are not accounted for. Included is the template relation obtained later in Section |3.4| (cf. Equation 
O). Hereafter we will use y = Mj — 5 log h and x — log W — 2.5. 



3. The Tully-Fisher Relation 

The most widely applied methods for the determination of redshift-independent distances of galaxies 
rely on the combination of photometric (distance dependent) and kinematic (distance independent) 
parameters. Among them, and arguably the most accurate, are the TF and fundamental plane relations, 
used for spiral and elliptical galaxies, respectively. With both these methods distances for individual 
galaxies, and therefore peculiar velocities, are computed by comparison with a fiducial template relation 
which must be observationally derived. The template relation defines the rest reference frame, against 
which peculiar velocities are to be measured. The importance of accurately calibrating such a tool cannot 
be overemphasized. It is possible that the discrepancies between different claims of bulk motions may be 
partly related to insufficiently well determined template relations. 

Ideally, the template employed is tied to a kinematical rest frame. In practice, the template originates 
from data within the sampled volume itself. We outline here our approach to calibrate the TF relation 
using our sample of data from 52 rich Abell clusters distributed throughout the sky. As we shall see, the 
TF zero point is obtained through an iterative process: the computation of the cluster incompleteness bias 
requires prior knowledge of the TF template zero-point, while the computation of the zero point demands a 
correction for the cluster incompleteness bias. The process does converge, however, in a small number of 
iterations. 
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3.1. The Calibration of the Tully-Fisher Relation 



For an assumed linear TF relation we need to determine two main parameters: a slope and a magnitude 
offset, or zero point. The slope of the TF relation is best determined by a sample that maximizes the 
dynamic range in (M,\ogW), i.e. one that preferentially includes nearby objects. On the other hand, the 
magnitude zero point of the relation is best obtained from a sample of distant objects for which a peculiar 
motion of given amplitude translates into a small magnitude shift. An often-adopted practice to calibrate 
the TF relation uses one relatively distant cluster of galaxies, for example the rich Coma cluster at cz~7200 
km s^ 1 . Such a choice relics on the assumption that galaxies within a cluster are essentially at the same 
redshift - any differences in their radial velocities are attributed to the cluster's virial stretch. Thus, they all 
equally participate in the local peculiar velocity field and they should all obey the same local TF relation. 
Another reason to choose a cluster like Coma for calibration is that, as mentioned above, the selection 
of a distant object limits the impact of the object's peculiar velocity V, at least to the extent that the 
physical size of the cluster is small in comparison to its mean distance; peculiar velocities introduce relative 
distortions of redshifts that are larger for nearby objects than for distant ones. In terms of magnitude, if 
V pcc <<cz, the zero point will be off by Am ~ — 2.l7V pcc /cz magnitudes (cf. Equation |l4| ). 

There are several problems with the above calibration scheme. First, a "template" cluster needs to have 
a large sample of spiral galaxies. Second, even for a cluster as distant as Coma, typical cosmic velocities 
may bias the relation's zero point. If the cluster were moving at a plausible speed of 500 km s _1 , then all 
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other estimates of peculiar velocities that use the template will systematically be off by 500(cz/czt em piate)> 
regardless of statistical uncertainties. 

To avoid such systematics in the TF relation, G97a,b use a "basket of clusters" approach. Their 
template is derived through an iterative procedure that simultaneously determines the TF zero point and 
the cluster motions of their sample. They assume that the mean peculiar velocity of the clusters farther 
than cz — 4000 km s _1 is null, i.e. they zero the monopole of the more distant cluster peculiar velocity 
distribution function. This approach does not, however, affect the value of the dipole or higher moments, 
and thus still allows an effective measure of possible bulk flows. The proximity of the SCI sample provides 
the stretch in galactic properties necessary to determine accurately the TF slope (it is for this reason 
that we adopt the slope computed in G97b). Also, the large database afforded by the SCI allows the 
statistical uncertainty of the magnitude zero point to be reduced to 0.02 — 0.03 magnitudes. The systematic 

1 /2 

uncertainty in the zero point is larger. For a sample of N objects with an rms velocity of (V 2 ) 1 at a mean 
redshift of (cz), the expected accuracy of the zero point is limited by systematic concerns to 

2.17(V^ 2 ) 1/2 

- - mag. (2) 



(cz) VN 

G97b conclude the SCI systematics only allows the universal TF zero point to be determined to within 0.04 
magnitudes. Thus the overall zero point accuracy obtainable from the SCI sample is ~0.05 magnitudes. We 
adopt the procedure described in G97b to calibrate the TF zero point using the SCII sample. The resulting 



improvement in its calibration is discussed in Section 3.4 



3.2. The Scatter of the Tully-Fisher Relation 

Any relation involving observed parameters has a limited accuracy described by the amplitude of the 
relation's scatter. Claims of the scatter in the TF relation vary from as low as 0.10 mag (Bernstein et al. 
1994) to as high as 0.7 mag (Sandage et al. 1994a, b; 1995; Marinoni et al. 1998). The amplitude of the 
scatter does depend on wavelength, and studies in the / band typically yield the tightest relations. The 
efforts with the largest samples yield la dispersion values of ~0.3-0.4 mag (Mathewson, Ford & Buckhorn 
1994; Willick et al. 1995; G97b). Uncertainties in observational measurements are not the only factors that 
lead to the overall spread in the data. The corrections to the observed fluxes and disk rotational velocities 
described in Papers I and II are not exactly known, nor are the methods we detail later to account for 
inherent sample biases such as cluster incompleteness. Moreover, there is an intrinsic component to the 
TF dispersion since individual galaxies have diverse formation histories. In fact, Eisenstein & Loeb (1996) 
advocate an intrinsic scatter of 0.3 magnitudes, a number greater than most estimates from observational 
work. In light of this fact, they make the interesting claim that either spirals formed quite early or that 
there must be a type of feedback loop that promotes galactic assimilation. 

As already established in G97b and Willick (1999), Figure ^reinforces the notion of low intrinsic 
scatter. The two dotted lines indicate the velocity width and magnitude uncertainties e x and e v , with 
e x multiplied by the TF slope b to put it on a magnitude scale; the solid line labeled e m = ^J(be x ) 2 + e y 
represents the average measurement uncertainty. The data displayed in Figure || are generated using equal 
numbers of galaxies per data point. The circles plotted represent the average standard deviations of the 
residuals from the fiducial TF relation. We see that the velocity width errors dominate those from the I 
band fluxes, which are approximately independent of velocity width. Furthermore, the logarithmic velocity 
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Fig. 5. — The observed scatter in the TF relation. Running averages (in units of magnitudes) of the 
observed errors are given as dotted lines. The lower solid line is the sum in quadrature of the two dotted 
lines and represents the total measurement uncertainty. The circles plotted indicate the standard deviations 
of residuals from the template TF relation. The upper solid line is the total scatter in the TF relation and 
has been computed by a sum in quadrature of the observed error distribution and an intrinsic component 
equivalent to that found in G97b (see Equation^ below). 

widths become increasingly uncertain for slower rotators (cf. G97b; Willick et al. 1997). We approximate 
the total observed scatter with a simple linear relation that depends on the velocity width: 

cr tot = -0.40a: + 0.38 mag. (3) 

The total scatter found here is in general larger than that found in G97b (crtot — —0.32a; + 0.32). This 
is unsurprising given our use of optical rotation curves instead of 21 cm profiles, a comparatively easier 
source from which to estimate velocity widths, and since the nearer SCI galaxies generally have better 
determined disk inclinations. The gap between the observed scatter and the measured errors is attributed 
to an intrinsic scatter contribution: the top line is a sum in quadrature of our observed measurement errors, 
e m , and the intrinsic scatter found in G97b: 

cr int = -0.28a; + 0.26 mag. (4) 

Simply put, our findings for an intrinsic component of the TF relation agree with those of Giovanelli and 
Willick and their collaborators. 



3.3. Observational Biases 

3.3.1. Cluster Population Incompleteness 

Quite possibly the most important selection effect to consider for this observing program is that of 
cluster incompleteness, the preferential sampling of the bright end of the luminosity function (LF). Several 



- 13 - 



authors have cautioned that cluster incompleteness can significantly alter inferred distances to clusters, 
though this ultimately depends on the amplitude of the scatter in the TF relation (see G97b, Willick 1994, 
Sandage, Tammann & Federspiel 1995, and the review by Teerikorpi 1997). Schechter (1980) first advanced 
the notion that the repercussions of this selection effect can be circumvented through an inverse fitting 
procedure, whereby the roles of the ordinate and abscissa in the TF diagram are reversed. The argument 
goes as follows: if there are no observational effects working against the selection of velocity widths and 
if the errors on the absolute magnitudes are negligible, then a fit to logVF vs. M will not be affected 
by a cutoff at faint magnitudes. Unfortunately, magnitude errors cannot be ignored (see Figure [5]) , and 
moreover, the errors in M and logVT" are coupled through inclination corrections; TF data do not obey a 
sharp faint magnitude limit. Inevitably, the absence of faint galaxies in a TF sample artificially brightens 
the zero point and lowers the slope. We describe next our methodology to account for the effect. 

Our first concern is to quantify the characteristics of the observed luminosity distribution with respect 
to the actual LF. As an aside, we note that current work by Andreon (1998) supports the notion of 
a canonical spiral LF by proposing that the LF for the separate morphological classes E, SO, and S is 
independent of environment - observed differences in the overall LF for different environments are merely 
due to varying proportions of the morphological classes (see, however, Iovino et al. 1999). Our first step is 
to sum up the observed luminosities in bins of absolute magnitude. Computing such absolute magnitude 
histograms for each individual cluster, where the membership counts can be as small as ~5, would not be 
statistically indicative on the whole. We therefore compute average histograms for several distance ranges of 
width 2000 km s _1 . We define the completeness function as the ratio of the observed luminosity distribution 
to the assumed intrinsic luminosity distribution. The following figures have been constructed assuming 
a Schechter LF with M* = -21.6 and a = -0.5 (see Figure 20 of Paper IV), but final TF templates 
using other Schechter parameters will also be provided. Figure ^ displays a smoothed representation of 
the completeness functions, along with our fitted approximation. As in G97b, we borrow a relation from 
Fermi-Dirac statistics to model completeness: 

C ^ = e (y-v!)/P + 1 (5) 

where yt represents a transition luminosity in the fit, and (3 characterizes the steepness of the drop. It should 
be noted that the final TF template is rather robust in terms of the completeness function construction. 
The template is largely insensitive to the choice of distance regimes, variations in the Fermi-Dirac fit profile, 
and the luminosity bin widths. Reasonable alterations in these parameters induce negligible changes in the 
TF template zero point. With the estimation of the completeness function in hand, we can now quantify a 
given cluster's incompleteness bias via Monte Carlo simulations. 



3.3.2. Monte Carlo Simulations of Cluster Parent Samples 

Our task is simplified by assuming the bias-corrected slope derived in G97b is valid for our sample, a 
justifiable assumption given the relative propinquity of their data - a broader dynamic range can be more 
easily sampled in nearby clusters. Furthermore, least-square fits to the data yield slopes that, within the 
errors, agree with the TF slope from G97b. On a more qualitative level, inspection of Figure ^ indicates 
our data agree with the slope from G97b. It is thus sufficient to only determine the bias in the TF offset 
for each cluster. 

We compute the bias as follows: for each galaxy a large number (iVit er = 10 3 ) of trial TF data points 
are generated that mimic the general characteristics of the actual data. Within each trial, a random 
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Fig. 6. — The distributions of cluster completeness for separate distance ranges. The dashed line is our 
smoothed approximation to the completeness function. The number of galaxies within each redshift range 
are indicated. 



(Gaussian deviate) magnitude offset Ay tr iai is first chosen according to the TF scatter relation, Equation 
|3|, which depends on the galaxy's measured velocity width parameter x. The addition of this offset to the 
magnitude inferred from the final, template TF relation yields the trial absolute magnitude: 

2/trial = hfX + d t f + Atrial- (6) 

The trial magnitude is kept only if it is a likely magnitude, i.e. a second random number from the interval 
[0,1] must be less than c(j/trial)> the completeness value for the trial luminosity. Otherwise, the process is 
repeated until a likely magnitude is found. After all Nn CI iterations are complete for a given galaxy, the 
incompleteness bias is taken to be the mean difference between the trial luminosities and that expected 
from the TF relation, i.e. 

Niter 

- Ay icb = ^ M ytriaM - (b tf x + otf). (7) 

■"'iter 

Figure M gives the biases calculated as a function of velocity width and the polynomials we fit to them. The 
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Fig. 7. — The incompleteness bias as a function of rotational velocity width is shown. Individual points are 
the biases generated for each galaxy from Monte Carlo simulations and the solid lines are polynomial fits to 
the bias. 

general property of the computed bias is as expected: incompleteness biases "turn on" at higher velocity 
widths for more distant objects. An incompleteness bias-corrected zero point for each cluster is then 
extracted from the cluster's distribution of observed velocity widths and bias-corrected absolute magnitudes. 
A tabulation of the offsets from the template zero point, a — a t { , is provided in Table |2| The cluster names 
are listed first and are followed by the number of cluster members N t { with reliable photometry and velocity 
widths. The next two columns are measures of the cluster offsets from the template zero point, the second 
of which is corrected for the effect of cluster incompleteness and includes an indication of its uncertainty e a 
in parentheses, e.g. —0.06(10) implies —0.06 ± 0.10 magnitudes. The last column of data listed is a a , the 
dispersion in the difference between the template and the cluster's set of absolute magnitudes. It is used to 
compute the offset uncertainty: 

e Q = er max /viVtf, (8) 

where er max is taken to be the maximum of [0.35 mag, a a ] to avoid overly optimistic measures of offset 
uncertainty in cases of chance alignment of the data due to small number statistics. Lastly, we remark that 
the proper incorporation of selection effect corrections is vital to determining the TF template, but the 
exact evaluation of such biases necessarily requires a TF template. A stable solution to this circular process 
is found within a few iterations. 
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Table 2: Cluster Offsets 
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N ti 


flbias — dtf 


a — dtf 


o a 


Cluster 


N 


fflbias — Otf 


a — dtf 


a a 


A2806 


10 


-0.15 


-0.13(11) 


0.27 


A3528 


3 


0.03 


0.18(07) 


0.12 


A 114 


9 


-0.01 


0.07(14) 


0.41 


A 1736 


6 


-0.03 


0.01(18) 


0.44 


A 119 


6 


0.04 


0.05(16) 


0.38 


A 1736b 


4 


-0.03 


-0.03(18) 


0.32 


A2877 


7 


-0.01 


0.03(15) 


0.40 


A3558 


8 


-0.13 


-0.10(15) 


0.43 


A2877b 


5 


-0.10 


-0.08(16) 


0.14 


A3566 


9 


-0.13 


-0.03(12) 


0.23 


A 160 


6 


-0.06 


-0.05(18) 


0.45 


A3581 


4 


-0.01 


0.04(20) 


0.40 


A 168 


9 


-0.13 


-0.12(13) 


0.38 


A 1983b 


8 


-0.27 


-0.26(12) 


0.25 


A 194 


13 


0.05 


0.09(13) 


0.45 


A1983 


7 


-0.14 


-0.07(19) 


0.50 


A 260 


9 


0.13 


0.23(14) 


0.42 


A2022 


8 


0.06 


0.14(12) 


0.32 


A 397 


14 


-0.17 


-0.13(15) 


0.57 


A2040 


10 


-0.10 


-0.03(14) 


0.43 


A3193 


6 


-0.14 


-0.10(14) 


0.09 


A2063 


18 


-0.18 


-0.14(09) 


0.36 


A3266 


2 


0.12 


0.31(25) 


0.21 


A2147 


19 


-0.09 


-0.06(09) 


0.39 


A 496 


9 


-0.18 


-0.13(12) 


0.36 


A2151 


22 


-0.08 


-0.06(08) 


0.40 


A3381 


4 


-0.17 


-0.16(18) 


0.12 


A2256 


8 


-0.04 


-0.01(12) 


0.33 


A3407 


8 


-0.07 


0.03(21) 


0.58 


A2295b 


4 


-0.05 


0.05(18) 


0.36 


A 569 


13 


0.02 


0.06(10) 


0.30 


A2295 


10 


0.07 


0.10(12) 


0.39 


A 634 


8 


0.04 


0.06(12) 


0.35 


A3656 


6 


0.01 


0.03(14) 


0.28 


A 671 


9 


-0.05 


0.02(12) 


0.34 


A3667 


4 


0.27 


0.37(18) 


0.28 


A 754 


3 


-0.07 


0.01(42) 


0.73 


A3716 


14 


-0.08 


-0.06(09) 


0.22 


A 779 


14 


0.01 


0.03(09) 


0.27 


A3744 


11 


-0.03 


0.03(11) 


0.37 


A 957 


6 


0.08 


0.13(14) 


0.16 


A2457 


9 


-0.02 


0.02(12) 


0.32 


A1139 


11 


-0.17 


-0.13(12) 


0.39 


A2572 


5 


-0.11 


-0.08(16) 


0.13 


Aim 


6 


-0.11 


-0.01(15) 


0.37 


A2589 


6 


-0.03 


0.04(14) 


0.23 


A1213 


6 


-0.12 


-0.12(14) 


0.33 
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12 


0.09 


0.13(10) 


0.30 


A1228 


13 


0.04 


0.12(10) 


0.33 


A2657 


5 


-0.03 


-0.01(16) 


0.26 


A1314 


8 


-0.04 


0.03(12) 


0.31 


A4038 


7 


-0.06 


-0.02(13) 


0.35 
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3.3.3. Homogeneous Malmquist Bias 

An observational bias commonly referred to as the "homogeneous Malmquist bias" is a distance 
underestimate for a galaxy drawn from a uniform distribution of galaxies. The bias is a direct result of 
the error on the measured distance - a galaxy taken from a Poissonian spatial distribution and measured 
to be at a distance modulus /i m ± A/i is more likely to actually lie between [i m + A^ than in the range 
yu m — A/i, due to the larger volume of the more distant shell. Following the reasoning given in Lynden-Bell 
et al. (1988), it can be shown that a measured distance R m is, on average, an underrepresentation of the 
true distance by the factor exp(3.5A 2 ): 

V = cz - H a R m ^cz- i? i? m e 3 ' 5A2 . (9) 

The factor A — 10 2ea — 1 is the fractional error in the TF distance measurement, roughly 17% for 
individual galaxies and 6% for our clusters. It was shown in G97b that such a correction for clusters of 
galaxies is rather small, almost to the point of being negligible (of order 1% on individual cluster distances). 
The inclusion of the above correction factor in later computations in this work does not change results 
appreciably. The influence of the homogeneous Malmquist bias on the inferred local bulk flow is discussed 
in Paper III. 

The Abell/ACO clusters do not represent a homogeneously distributed population. Thus, the 
application of a homogeneous Malmquist bias correction may appear incorrect. However, a correction 
that takes into account the clustering properties of clusters differs from the homogeneous Malmquist bias 
correction only in the second order. Given the miniscule size of the correction and the uncertainty of 
a possible "inhomogeneous" Malmquist bias correction, we consider the application of a homogeneous 
Malmquist bias correction a satisfactory approach. 

3.3.4. Cluster Size 

In calculating absolute magnitudes, we assume that each cluster member within one Abell radius of 
the cluster center is at a distance equal to the average of all the redshifts available for the cluster from the 
literature. It is unlikely that a given galaxy will actually be at this distance, but it is a useful approximation 
if the distance to the cluster is significantly larger than the cluster's virial size. For all but a handful of our 
clusters this is the case, so cluster size biases play a small role in our analysis. 

Our concern here is the method of averaging used. The calculation of a cluster offset involves averaging 
over absolute magnitudes. Thus, averages of the logarithms of the distances are computed, when the 
logarithm of the average distance is actually desired. This results in a systematic underestimate of cluster 
distances. We can use either the angular distribution on the sky of the cluster galaxies or the Abell radius 
to infer the approximate physical size of the each cluster. Simple analytic calculations using either estimate 
of the cluster size show that the amplitude of the bias is at most of order 0.001 magnitudes for even our 
closest clusters, so we shall not concern ourselves with this effect any further. A more serious concern, the 
morphological dependence on the TF relation, is investigated next. 
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3.3.5. Morphology 

Ideally, scientific models are simple and are thus constructed with the fewest number of free parameters 
the data demand. In TF work, it would be preferable to limit the range of morphologies sampled, since 
physical parameters are known to vary along the Hubble sequence (Roberts & Haynes 1994). For example, 
it has been shown that, for a given (optical) luminosity distribution, the velocity width distribution for 
early-type galaxies is shifted to higher rotational speeds with respect to the distribution for later types 
(Roberts 1978). There is evidence, however, that the TF differences between morphological types appear to 
diminish at wavelengths longer than / band (e.g. Aaronson & Mould 1983). Regardless, we are fortunate 
that one sample selection criterion of ours, that of disks being rich in ionized gas, encouraged a rather 
homogeneous sample comprised mainly of Sc types. Our sample has the following population properties: 
14% are of type earlier than Sb, 16% are type Sb, and 70% are classified as Sbc or later. The majority 
(52%) of the galaxies are type Sc. 

Our sample affords statistically significant tests of TF morphological dependencies. We plot in Figure 
H the TF parameters with symbols differentiated according to morphological class. Filled circles symbolize 




2.1 2.4 2.7 



Fig. 8. — The TF diagram differentiated according to morphological class. Asterisks represent galaxies 
earlier than Sb, filled circles symbolize type Sb, and open circles are for Sbc and later. The data are 
corrected for cluster peculiar motions. 



Sb types, whereas open circles and asterisks represent types later and earlier than Sb, respectively. The 
plot includes the data for all galaxies deemed to have reliable velocity widths, with each cluster member 
corrected for cluster peculiar motion. The solid line drawn uses the fiducial TF slope of 6 t f = —7.68 and 
zero point a t f = —20.91 mag (cf. Section 3.4). As found in G97b, a clear distinction is evident between the 
three Hubble types in the form of a fainter zero point for earlier types. The error-weighted averages of the 
offsets Amy from the template zero point for the three morphological classes differ in the following ways 
for our sample: 



Types earlier than Sb: —0.27 (—0.32) mag 
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Type Sb: -0.11 (-0.10) mag 
Other types: unchanged 

The numbers given in parentheses are those from G97b. For consistency, we will continue to utilize the 
offsets obtained in G97b. 



3.3.6. Environment 

Yet another possible bias in our sample is the effect of environment. For instance, the more distant 
clusters in our sample (and the Abell/ACO catalog in general) tend to be richer. A richer cluster typically 
has a stronger intracluster medium X-ray flux and a more regular, elliptical-dominated core (see, for 
example, Sarazin 1986). Spiral galaxies predominantly lie in a rich cluster's periphery and the closer a spiral 
disk is to the cluster center, the less likely it is to contain neutral hydrogen gas (Giovanelli & Haynes 1985). 
This lack of interstellar gas within clusters of galaxies may be due to evaporation into the hotter intracluster 
gas, or it may be attributed to "stripping" originating from either tidal galaxy-galaxy interactions or ram 
pressure ablation on intracluster gas. Rubin, Ford, & Whitmore (1988) and Whitmore, Forbes, & Rubin 
(1988; WFR hereafter) claim that the inner spiral galaxies within clusters exhibit falling rotation curves 
(RCs), as opposed to the asymptotically flat or rising RCs usually seen in the cluster periphery and field. 
Furthermore, they find cluster RCs may be of lower amplitude than field RCs. They offer the explanation 
that the falling (and lower amplitude) RCs arise from mass loss — the inner galaxies have had their dark 
matter halos stripped — or that the cluster environment simply inhibits halo formation. A related finding 
by WFR is a monotonic increase in the mass to light ratio, with cluster radius, which they ascribe to the 
changing RC shape with cluster radius. This view has been contested, however, by Amram et al. (1993) 
and Vogt (1995) who find little evidence for any gradients in the outer portions of RCs. If RCs in the inner 
regions of rich clusters do indeed differ from RCs in other environments, the ramifications are significant. 
The dependence of the TF relation and/or its dispersion on environment are possible consequences. Finally, 
we point out that rich clusters represent high density peaks and may thus be home to a recent merger; 
at least one third of all rich clusters are home to a recent inhomogeneous superposition of two or more 
separate systems (Girardi et al. 1997). It is prudent to verify whether environment plays a significant role 
in our data. 

Figure || displays one such test of environmental bias. For each cluster member we have plotted 
residuals from the TF relation as a function of projected distance from the nominal cluster center. Panel 
(a) displays the data for all cluster members, while panels (b) and (c) differentiate between Abell cluster 
richness. All data are corrected for peculiar velocity and morphological offsets. We see no apparent trend 
with projected cluster radius in any of the panels. This is in agreement with the work of Biviano et al. 
(1990) and G97b, which showed no change in the TF relation for different environments. We consequently 
do not consider environmental bias to be a serious concern with our data. 



3.4. The Template Relation 

Let afc be the TF zero point for cluster k and Ay pcc k = ak — atf be the shift due to the cluster's 
peculiar motion. The i th cluster member's absolute magnitude can then be expressed as 



Ui — Vcor,i 



(10) 
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"(a) All 



(b) R = 25 clusters: 




"(c) R>0 27 clusters 



2 - • .it— 



-1 



12 3 

hr (Mpc) 

Fig. 9. — A test of environmental bias is displayed. TF residuals are given versus projected physical distance 
from the cluster center. The data are corrected for peculiar velocity and morphological offsets. Panel (a) 
displays all the cluster data, whereas panels (b) and (c) differentiate between cluster richness. 



where Aj/i c b is the cluster's incompleteness bias correction described in Section 3.3.1 and 

51og 10 (cz c ius/100) — 25, TOcor being the corrected apparent magnitude given by 
Equation 1 in Paper II and cz c i us is the cluster systemic recessional velocity in units of km s . We compute 
the fc th cluster's zero point by averaging over the individual cluster members, i.e. 



(ex,i6) ! 



(11) 



with Nk representing the number of cluster members in the k th cluster.^] The uncertainty e a ^ of the above 
zero point computation is described by Equation |[ The calculation of the template zero point Otf follows. 
Alternative estimators of att are: 



Otf 



Ek^k 



EfiA 



(12) 



The first two estimators are averages over the individual cluster zero points, weighted either by the number 
of cluster members or by the cluster zero point uncertainties. The third calculation is a simple average over 
all N cluster galaxies, with each galaxy weighted by its total uncertainty. The three estimators give global 
zero points that agree within the estimated errors (0.02 magnitudes); we adopt the third computation. Our 



5 Though we are only solving for the relation's offset here, this mode of calculation is similar in spirit to the bivariatc 
calculations described in G97b, in that both magnitude and velocity width errors are considered. 
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definition of a TF template zero point is based on the assumption that the average peculiar velocity of the 
cluster set is null. 

Table || lists results for several subsets of the data and for two different Shechter LFs: 
a = -0.50, M* = -21.6 and a = — 0.75, M* = —22.0. A few fit parameters are listed, including 
the number of galaxies used, the TF template zero point and its associated total uncertainty (i.e. 
considering both statistical concerns and the kinematic uncertainty described by Equation ^) , an estimate 
of the scatter, and chi squared divided by the number of degrees of freedom. Notice that the dispersion 
and the zero point uncertainty remain relatively unchanged across all subsets. The fact that the zero 
point differs slightly for late and early type galaxies is expected - we saw earlier that the morphological 
type correction of —0.32 magnitudes advocated in G97b for Sa and Sab galaxies is possibly too large for 
our sample by a few hundredths of a magnitude. We should also point out that the slightly brighter zero 
point seen for objects beyond 10,000 km s _1 works against recent claims by Tammann (1998) and Zehavi 
(1998) that the Hubble constant is higher within ^10,000 km s _1 than it is beyond this distance. For our 
purposes, a fractional decrease by cr in H a beyond 10,000 km s _1 should yield a zero point that is fainter 
by 6h5 log 10 e, or 0.04 mags for en = 0.02. In contrast, our data reflect a slightly brighter zero point beyond 
10,000 km s _1 . 

Figure [Io| gives the TF plots for each cluster and Figure |ll| combines the data for all SCII galaxies. 
The data are corrected for cluster incompleteness bias and cluster peculiar motion, in accordance with 
Equation |l0|. In the A2877, A1736, A1983, and A2295 panels, the error bars containing filled circles 
represent members of "A2877b," "A1736b," "A1983b," and "A2295b," respectively. The new TF template 
is drawn as well: 

Mi — 5 log/i = — 7.68(log W — 2.5) — 20.91 mag. (13) 



The residuals given in Figure [L2| indicate the quality of the cluster membership assignments. Since the 
abscissa is the difference between individual and cluster redshifts, the center of each cluster corresponds to 
the point (0,0). Moreover, residuals from field galaxies that are receding according to Hubble expansion but 
were incorrectly declared cluster members should, on average, lie on the line of slope 5. There is no clear 
evidence for a significantly large subsample of improperly assigned memberships. 



4. The Peculiar Velocity Sample 

We interpret the departure of a cluster's average zero point from that of the template relation as an 
indication of peculiar motion, with larger departures from the template implying larger amplitude peculiar 
velocities. Quantitatively, for a cluster at a redshift z with an average departure from the template of 
a — <Ztf, we write the peculiar velocity as 

V = az(l-10 - 2(o_atf) ). (14) 

We compute peculiar velocities and a measure of their errors for the SCII cluster sample and display the 
results in Table || and Figure [l^, an Aitoff projection of Galactic coordinates. The symbols plotted in the 
figure reflect both the radial directions of the peculiar velocities and the strengths of the measurements - in 
the CMB reference frame, open (filled) circles represent approaching (receding) clusters and the circle size 
is inversely proportional to the accuracy of the measurement. The largest cluster peculiar velocities, e.g. 
those for A3266 and A3667, are also the most uncertain, as the clusters are poorly sampled. 
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Table 3: Template Fit Parameters 



Sample 


N 


a 






x 2 / dof 


no incompleteness bias 


441 


-20.96 


0.02 


0.39 


1.1 


qlp = -0.50, M* = -21.6 














441 


-20.91 


0.02 


0.38 


1.1 


2.5cr clip 


425 


-20.91 


0.02 


0.34 


0.9 


60° < i < 80° 


252 


-20.91 


0.02 


0.38 


1.2 


i > 45° 


429 


-20.91 


0.02 


0.38 


1.1 


W > 150, y < -18.5 


428 


-20.92 


0.02 


0.37 


1.1 


T < Sb 


139 


-20.95 


0.03 


0.38 


1.3 


i > SO 


302 


-20.89 


0.02 


0.37 


1.0 


cz > 10, 000 km s" 1 


325 


-20.92 


0.02 


0.38 


1.1 


-^mcm ^ 5 


413 


-20.92 


0.02 


0.38 


1.1 


^proj <1 Abell radius 


347 


-20.91 


0.02 


0.39 


1.1 


\cz — cz c i us | < 1000 km s _1 


368 


-20.91 


0.02 


0.38 


1.1 


includ. foreground & background 


510 


-20.90 


0.02 


0.39 


1.2 


qlp = —0.75, M* = -22.0 














441 


-20.90 


0.02 


0.38 


1.1 


2.5cr clip 


425 


-20.90 


0.02 


0.34 


0.9 


60° < i < 80° 


252 


-20.90 


0.02 


0.38 


1.2 


i > 45° 


429 


-20.90 


0.02 


0.37 


1.1 


1¥ > 150, y < -18.5 


427 


-20.91 


0.02 


0.37 


1.1 


T < Sb 


139 


-20.95 


0.03 


0.38 


1.3 


T> Sb 


302 


-20.88 


0.02 


0.37 


0.9 


cz > 10, 000 km s" 1 


325 


-20.91 


0.02 


0.38 


1.1 


^Vmem ^ 5 


413 


-20.91 


0.02 


0.34 


1.1 


^proj <1 Abell radius 


347 


-20.90 


0.02 


0.38 


1.1 


cz — C2; c i us | < 1000 km s _1 


368 


-20.90 


0.02 


0.38 


1.1 


includ. foreground & background 


510 


-20.89 


0.02 


0.39 


1.2 
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Fig. 10. — TF plots for the clusters are given. Only cluster members are included in these plots. Each 
luminosity is corrected for cluster incompleteness and cluster peculiar velocity. The morphological offset 
described in Section 3.3.5| for Sb and earlier type galaxies is incorporated. The new TF template relation is 



drawn as a reference (Equation 13) 



Our sample includes the central portions of various high density peaks and/or superclusters in the 
local Universe. It is significant to note that their motions are consistent with small departures from rest in 
the CMB frame, within the quoted error. For instance, the Shapley Supercluster, represented here by its 
core (A3558) and some peripheral members (A1736, A3528, and A3566), has an average CMB velocity of 
118±495 km s _1 . The Hercules region (A2147 and A2151) and the A2572/2589/2593/2657 supercluster are 
also slow movers, with average peculiar velocities of 307 ± 301 km s _1 and —222 ± 372 km s _1 , respectively. 
These small motions are consistent with the notion that such massive systems best represent "kinematic 
anchors" in the local velocity field. 
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Figure HO (Continued) 
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4.1. The One-Dimensional Peculiar Velocity Distribution 

It is useful to estimate the line-of-sight distribution of peculiar velocities. The amplitude of that 
distribution has been known to be a very sensitive discriminator of cosmological models (see, for example, 
Bahcall & Oh 1996). The SCII cluster sample is relatively distant and the cluster membership counts are 
relatively anemic when compared to the SCI sample of G97a. Consequently, SCII peculiar velocities are 
much less certain and the overall distribution shown in Figure [l4| is significantly broadened by measurement 
errors. The peculiar velocities are represented by equal area Gaussians centered at the peculiar velocity 
of each cluster with dispersions equal to the estimated peculiar velocity errors. The thick dashed line 
superimposed on the plot is the sum of the individual Gaussians (its amplitude has been rescaled for 
plotting purposes). The la dispersion in the observed distribution of peculiar velocities is found from a 
Gaussian fit to the dashed line: cid.obs = 796 km s _1 . This value, however, is biased high by measurement 
errors. Recovering an estimate of the true value can easily be obtained via Monte Carlo simulations, 
yielding aid — 341 ± 93 km s _1 where the error estimate derives from the scatter in the dispersions of the 
simulated samples. This value of aid is consistent with a relatively low density Universe (Giovanelli et al. 
1998b; Bahcall & Oh 1996; Borgani et al. 1997; Watkins 1998; Bahcall, Gramann & Cen 1994a; Croft & 



2.2 



2.4 2.6 



Fig. 11. — Plotted are the TF data for all galaxies, corrected for peculiar motions, incompleteness bias, 



and the morphological type offset described in Section 3.3.5. The new TF template relation is drawn as a 
reference (Equation Il3|) . 



Efstathiou 1994). 



5. Summary 

We have presented TF data and estimated the peculiar velocities of 52 rich Abell clusters spread across 
the sky and distributed between ~50 and 200h~ 1 Mpc. Optical rotation curves and I band photometry for 
522 spiral galaxies in the fields of these systems have been obtained and presented in separate publications. 

In conjunction with the robust TF slope extracted from the relatively nearby SCI cluster sample of 
Giovanelli and coworkers, we find the / band TF relation to follow 

Mi - 5 log 10 h = -7.68(log 10 W - 2.5) - 20.91 mag. (15) 

The relation has an average scatter of 0.38 magnitudes. The zero point of the TF template has a statistical 
accuracy of 0.02 mag; combined with a kinematical uncertainty of 0.01 mag, which is limited by the 
assumption that the 52 clusters' average peculiar velocity is null, the overall uncertainty of the TF zero 
point is 0.02 mag. 

Peculiar velocities are obtained for each of the 52 clusters, with reference to the global template relation. 
The typical uncertainty on the peculiar velocity of each cluster is ±0.06cz, where cz is the mean cluster 



-0.1 0.1 

log 10 V - log 10 V(cluster) 

Fig. 12. — A test of cluster membership is displayed in the form of TF residuals as a function of separation 
from the cluster redshift. A line of slope 5 is given, the relation on which field galaxies incorrectly assumed 
to be cluster members should lie. 




Fig. 13. The all-sky distribution of the SCII peculiar velocity sample in Galactic coordinates. The symbol 
diameters are inversely proportional to the peculiar velocity uncertainties. The two examples in the lower 
left give the scale. Filled (open) circles refer to positive (negative) peculiar velocities in the CMB frame. 
Asterisks mark the apices of the motion of the Local Group with respect to the CMB and the LP cluster 
inertial frame. 

velocity. The rms line-of-sight component of the cluster peculiar velocity for our sample, dcbroadened for 
measurement errors, is 341 ± 93 km s _1 . This number agrees with that determined for the SCI sample. 
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Fig. 14. — A display of equal area Gaussians representing the peculiar velocity sample. Each Gaussian is 
centered at the value of a cluster peculiar velocity, has a dispersion given by the uncertainty in the peculiar 
velocity measurement, and has an amplitude that is inversely proportional to the uncertainty. The sum of 
the Gaussian profiles is given by the thick dashed line. 

The results presented here are based on observations carried out at the Palomar Observatory (PO), at 
the Kitt Peak National Observatory (KPNO), at the Cerro Tololo Inter- American Observatory (CTIO), 
and the Arecibo Observatory, which is part of the National Astronomy and Ionosphere Center (NAIC). 
KPNO and CTIO are operated by Associated Universities for Research in Astronomy and NAIC is 
operated by Cornell University, all under cooperative agreements with the National Science Foundation. 
The Hale telescope at the PO is operated by the California Institute of Technology under a cooperative 
agreement with Cornell University and the Jet Propulsion Laboratory. This research was supported by 
NSF grants AST94-20505 and AST96-17069 to RG and AST95-28960 to MH. LEC was partially supported 
by FONDECYT grant #1970735. 
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Table 4: Cluster Peculiar Velocities 



Cluster 


a 


- Otf 


CZhol 


C2i g 




Vtf.cmb 


^peCjCmb 


A2806 


-0.135 


112) 


8019 


7918 


7867 


7403 


464 


0382) 


A114 


0.071 


134) 


17436 


17487 


17144 


17722 


-578 


1111) 


A119 


0.053 


153) 


13470 


13607 


13141 


13416 


-275 


;0988) 


A2877 


0.039 


155) 


7165 


7097 


6974 


7078 


-104 


;0489) 


A2877b 


-0.030 


158) 


9231 


9163 


9040 


8733 


307 


0634) 


A160 


-0.042 


176) 


12390 


12577 


12072 


11792 


280 


0977) 


A168 


-0.106 


127) 


13365 


13496 


13049 


12370 


679 


0725) 


A194 


0.089 


127) 


5343 


5460 


5037 


5253 


-216 


;0302) 


A260 


0.218 


162) 


10925 


11135 


10664 


11839 


-1175 


0835) 


A397 


-0.144 


149) 


9803 


9916 


9594 


9041 


553 


0630) 


A3193 


-0.100 


144) 


10558 


10381 


10522 


10072 


450 


0668) 


A3266 


0.392 


248) 


17774 


17564 


17782 


20482 


-2700 


2345) 


A496 


-0.128 


123) 


9860 


9779 


9809 


9243 


566 


0513) 


A3381 


-0.157 


176) 


11410 


11196 


11510 


10712 


798 


0868) 


A3407 


0.090 


221) 


12714 


12450 


12861 


13040 


-179 


1235) 


A569 


0.061 


099) 


5927 


5998 


6011 


6168 


-157 


;0280) 


A634 


0.088 


125) 


7822 


7917 


7922 


8144 


-222 


;0469) 


A671 


0.124 


118) 


15091 


15046 


15307 


15427 


-120 


;0838) 


A754 


0.047 


414) 


16282 


16052 


16599 


16691 


-92 


;3294) 


A779 


0.026 


095) 


6967 


6929 


7211 


7311 


-100 


;0320) 


A957 


0.119 


144) 


13464 


13261 


13819 


14685 


-866 


;0974) 


A1139 


-0.107 


118) 


11851 


11665 


12216 


11522 


694 


0629) 


A1177 


-0.008 


150) 


9755 


9661 


10079 
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